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Abstract. Variation in the adhesive protein gene se¬ 
quences of Mytilus edulis, M, galloprovincialis, and Al. 
trossulus collected in Delaware, Kamaishi (Japan), and 
Alaska, respectively, was analyzed by the polymerase 
chain reaction (PCR) using two sets of oligonucleotide 
primers. The first set. Me 13 and Me 14, was designed to 
amplify the repetitive region. The length of the amplified 
fragments was highly variable, even among samples of the 
same species. Another set. Me 15 and Me 16, was designed 
to amplify a part of the nonrepetitive region. The length 
of the amplified fragments was uniform in each species 
and differed interspecifically; 180, 168, and 126 bp for Al. 
edulis , Al. trossulus, and Al galloprovincialis, respectively. 
The amplified sequence of Al trossulus resembled that of 
Al edulis. Mussels from other sites were also examined 
by PCR using Me 15 and Me 16. Wild mussels from 
Tromso (Norway) and cultured mussels from Brittany 
(France) were identified as Al. edulis. Cultured mussels 
from the Mediterranean coast of France and wild mussels 
from Shimizu (Japan) were identified as Al. galloprovin¬ 
cialis. Some wild mussels from Hiura (Japan) were iden¬ 
tified as a hybrid between Al. galloprovincialis and Al. 
trossulus. Thus, the length of this part (variable region) 
of the sequence is proposed as a diagnostic marker for 
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these three morphologically similar species and their 
hybrids. 

Introduction 

Two types of polyphenolic proteins, foot proteins I and 
2. both of which incorporate 3,4-dihydroxyphenvlalanine 
(DOPA) into their primary structures, have been isolated 
from the mussel Mytilus edulis, and characterized (see 
Waite, 1992, for a review). Foot protein 1 is an adhesive 
protein that contains repeats of the decapeptide motif 
AKPSYP*P*TY*K. where P* and Y* denote hydroxy- 
proline and DOPA, respectively, and the hexapeptide 
motif AKPTY*K (Waite and Tanzer, 1981; Waite, 1983; 
Waite cl al.. 1985). Foot protein 2 is an epidermal growth- 
factor-like protein that forms the adhesive plaque matrix 
(Rzepecki el al.. 1992; Inoue el al.. 1995a). The primary 
structure of the foot protein 1 has been determined by 
cDNA or gene cloning. It has been shown that it consists 
of a relatively short nonrepetitive domain and a long re¬ 
petitive domain. The repetitive domain contains more 
than 70 decapeptide repeats and 13 or 14 hexapeptide 
repeats, but the number and distribution pattern of the 
motifs arc variable even in the same species (Filpula el 
al.. 1990; Laursen, 1992). 

We isolated cDNA encoding foot protein 1 from M. 
galloprovincialis sampled in Japan and compared the 
predicted amino acid sequence to those of M. edulis (Inoue 
and Odo, 1994; Inoue el al.. 1995b). The M. galhpro- 
vincialis sequence also consisted of the nonrepetitive and 
repetitive domains, but remarkable differences were ob¬ 
served in both. The major difference in the repetitive do- 
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Figure 1. Positions of PCR primers in the adhesive protein genes of Mvtilis edulis and M galloprovincialis. 
Adhesive protein genes are shown as cDNA sequences. Position and direction of primers is indicated by 
arrows. Solid and open boxes indicate the signal peptide and the repetitive region, respectively. Dotted boxes 
indicate the nonrcpetitive region, and the box with dense dots indicates the 18 amino acid sequence found 
in the M edulis sequence but not in that of M galloprovincialis. Solid lines indicate untranslated region. 
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main is that the M galloprovincialis sequence contained 
62 decapeptide repeats but no hexapeptide motif. In con¬ 
trast, the nonrepetitive domains differ by a deletion of 18 
amino acids observed in the M. galloprovincialis sequence. 

In this study, we prepared two sets of oligonucleotide 
primers for polymerase chain reaction (PCR) to amplify 
the whole repetitive region and a part of the nonrepetitive 
region. Using these, we have analyzed foot protein 1 se¬ 
quences of M. edulis , M. galloprovincialis , and M. tros- 
sulus. These three species are closely related (Gosling, 
1984; Koehn, 1991; Gardner, 1992; Seed, 1992, for re¬ 
views) and practically indistinguishable by morphological 
characteristics only. We report that the length of the frag¬ 
ments amplified from the nonrepetitive region is specific 
to each species, but the length of the repetitive region is 
highly variable even within the same species. Thus the 
nonrepetitive region can be used as a diagnostic marker 
for identification of the three species. 

Materials and Methods 

Mussels 

M edulis was collected at Lewes (Delaware, USA); M. 
galloprovincialis and M. trossulus were sampled at Ka- 
maishi (lwate, Japan) and Juneau (Alaska, USA). All these 
sampling points are “pure sites" of each species where 
the other two species do not occur (McDonald et al., 
1991), Wild mussels were also collected at Tromso (Nor¬ 
way), Hiura (Hokkaido, Japan), and Shimizu (Shizuoka, 
Japan). Mussels cultured in Brittany (France) and on the 
Mediterranean coast of France were obtained at a fish 
market at Ferncy-Valtaire (France). 

DNA extraction 

A piece of the gill from each mussel, about 0.5 cm 2 , 
was incubated in 500 jul lysis buffer containing 50 m/1/ 


Tris-HCl (pH 7.5), 10 m M EDTA, 0.5% SDS, 500 pg/m\ 
Proteinase K at 55°C for 2-4 h. Samples were then ex¬ 
tracted twice with equal volumes of saturated phenol and 
twice with phenol chloroform; isoamyl alcohol (24;24:1). 
The aqueous phase was precipitated with ethanol and dis¬ 
solved in 50 ^1 TE (10 m M Tris, 1 m M EDTA). 

PCR amplification 

About 100 ng of DNA was dissolved in 100jul 1 X Tth 
buffer (TOYOBO, Japan) containing 6 jug sense primer, 
6 pg antisense primer, and 200 pM dNTP. After pre¬ 
heating to 95°C, 1 unit of Tth DNA polymerase (TO¬ 
YOBO, Japan) was added and 30 cycles of amplification 
were performed. Each cycle consisted of 30 s at 94°C, 30 
s at 56°C, and 90 s at 70°C. The sequences of the primers 
were Me 13, CCA CTT GCA AAG AAG CTG TCA TCT; 
Me 14, ACA A AC GTT AAA ATG TGT AGT ACA 
GTA; Me 15, CCA GTA TAC AAA CCT GTG AAG A; 
Me 16, TGT TGT CTT A AT AGG TTT GTA AGA. 
Positions of primers in the foot protein 1 cDNA sequence 
in \L galloprovincialis are shown in Figure 1. 

Electrophoresis of amplified products 

Ten microliters of PCR product was mixed with the 
loading dye solution containing bromophenol blue (BPB) 
and xylene cyanol and subjected to agarose gel electro¬ 
phoresis. For analysis of PCR products, 4.8% NuSieve 3; 
1 agarose (FMC) was used for the nonrepetitive region, 
and 1% LE agarose (FMC) was used for the repetitive 
region. Electrophoresis on a 4.8% gel was continued until 
BPB reached the end of the gel. 

Sequencing 

The fragment amplified from the genome of M. tros¬ 
sulus using primers xMe 15 and Me 16 was isolated and 
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Figure 2. Representative results of amplification of the repetitive 
region of the adhesive protein gene. Amplified products were electro- 
phoresed on 1% agarose gel. Lanes l and 2, Mvtilis edulis; lanes 3 and 
4, M. irossulus; lanes 5 and 6, M galloprovincialis. M, molecular marker 
(lambda DNA digested with EcoT 141). 


inserted into the Snui I site of pUC19. Sequences of both 
strands of three independent clones were determined using 
a 373A DNA sequencer (Applied Biosystems Inc.) and a 
PRISM Dyeterminator Cycle Sequencing Kit (Applied 
Biosystems Inc.). 


Results 

)'aviation in the repetitive region 

The primers Me 13 and Me 14 were designed to amplify 
the repetitive region using the sequences identical to both 
M . edulis and M. galloprovincialis. Since the sense primer. 
Me 13, corresponds to a part of the nonrepetitive region 
and the antisense primer. Me 14, to a part of the 3' un¬ 
translated region, the whole repetitive region is amplified 
by PCR. M. edulis, M, galloprovincialis, and M. trossulus 
were collected at Delaware, Kamaishi, and Juneau, re¬ 
spectively. These sampling points are known to be “pure 
sites” at which no other species of the M. edulis complex 
is found (McDonald et ai, 1991). We analyzed 8, 16, and 
8 individuals of M. edulis, M galloprovincialis, and M . 
trossulus, respectively, using primers Me 13 and Me 14. 
Since the repetitive region is relatively long and highly 
repetitive, it was difficult to amplify the whole repetitive 
region if the template DNA was insufficiently pure and 
long, but prominent bands were successfully obtained by 
using well-purified, high molecular weight DNA. Typical 
results are shown in Figure 2. Sizes of the band ranged 
from 2.2 to 2.8 kb. The fragments obtained from M. edulis 
were generally" larger than those of the other two species. 
The sizes of bands in M. trossulus and M. galloprovincialis 
were similar but, on average, the former were slightly 
larger. Many individuals had two-banded (heterozygous) 


patterns, as expected for a highly variable polymorphism. 
One sample of M galloprovincialis exhibited three bands, 
which may be a naturally occurring triploid or a mosaic 
individual that possesses a cell lineage having the differed 
length of foot protein 1 gene. It is, however, also possible 
that the third band is a heteroduplex of two different frag¬ 
ments. 

1 filiation in the nonrepetitive region 

Another set of primers. Me 15 and Me 16, was also 
prepared to amplify a part of the nonrepetitive region 
using sequences perfectly identical between M. edulis and 
AI. galloprovincialis (Fig. 1). The size of the amplified 
fragment estimated from sequence data previously re¬ 
ported (Filpula et ai, 1990) is 180 bp in AI. edulis. In AI. 
galloprovincialis, the expected size is 126 bp because the 
sequence of AI. galloprovincialis contains a deletion of 18 
amino acids (Fig. 1; see also Jnoue and Odo, 1994). Using 
these primers. 8. 32, and 16 individuals of AI. edulis, AI. 
galloprovincialis, and AI. trossulus —including the same 
samples used in the analysis of the nonrepetitive region— 
were examined. PCR analysis indicated that all samples 
exhibited a single band. Representative results are shown 
in Figure 3. The position of the band was uniform in each 
species but differed from species to species. The size of 
the amplified fragments of AI. edulis and AI. gallopro¬ 
vincialis estimated by mobility in agarose agreed with 
those expected. Fragments from AI. trossulus were shorter 
than those of AI. edulis but longer than those of M. gallo¬ 
provincialis. To determine the length and sequence of the 
amplified fragment of AI. trossulus, the band obtained 
from one sample (Fig. 3, Lane 3) was isolated and sc- 



Kigurc 3. Representative results of amplification of lhe nonrepetitive 
region of the adhesive protein gene. Amplified products were electro- 
phoresed on 4.8% NuSieve GTG agarose gel (FMC). Lanes l and 2. 
Mvtilis edulis; lanes 3 and 4. At trossulus; lanes 5 and 6. M gallopro- 
vineiahs. M, molecular marker (pUC19 DNA digested with //c//>Il). 
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C CAGTATACAAACCTGT G AAGACAAGTTATT CGTCACCATATAAACCACCAACATACCAA 
PVYKPVKTSYSSPYKPPTYQ 

70 80 90 100 110 120 

CCACT CAAAAAGAAACCGATGGACTATAATAGTTCTCCGCCAACATATGGATCAAAGACA 
PLKKKPMDYNSSPPTYGSKT 

130 140 150 160 

AACT AT CTTGCAAAGAAGCTGTCA TCTTACAAACCTATTAAGACAACA 
NYLAKKLSSYKPIKTT 

Figure 4. Nucleotide and deduced amino acid sequences of the Fragment amplified from the genomic 
DNA of Mytilis irossulus using primers Mel 5 and Me 16. Underlined sequences were derived from primers. 


qucnced. The sequenced fragment including the primer 
sequences was 168 bp, 12 bp shorter than the correspond¬ 
ing region of M. cdulis (Fig. 4). This difference is small 
but measurable by mini-electrophoresis, as shown in Fig¬ 
ure 3. The nucleotide sequence of the amplified fragment 
was compared with corresponding sequences of M. cdulis 
and M. galloproviucialis (Fig. 5). Since the region shown 
in Figure 5 was especially variable among the three species, 
this region is hereafter referred to as the “variable region.” 
It seems that variation among the three species was caused 
by deletion or addition of short sequences, not by base 
substitutions. 

The variable region of mussels cultured in Brittany, 
those cultured on the Mediterranean coast of France, 
and the wild mussels collected at Tromso and Hiura 
were also examined by PCR using Me 13 and Me 16. 


Four individuals were examined in each group and rep¬ 
resentative results arc shown in Figure 6. All the wild 
Tromso mussels and all the cultured Brittany mussels 
exhibited the 180-bp fragment and were identified as 
M cdulis. All the cultured Mediterranean mussels and 
the wild mussels from Shimizu exhibited a 126-bp frag¬ 
ment and were identified as M. galloprovincialis. These 
results are consistent with the distribution map of mus¬ 
sels that was made by using allozyme characters 
(McDonald cl al. t 1991). Eight wild mussels of Hiura 
were also examined using Me 15 and Me 16. Two bands, 
126 and 168 bp, were amplified from six individuals 
(Fig. 6), but only the 126-bp band was amplified from 
the remainder (data not shown). Hiura is on Hokkaido 
Island, where both native M. irossulus and introduced 
M. galloprovincialis are distributed. The mussels that 
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Me 

Mg 


CAAGTT ATT CGT CACCAT AT AAACCACCAACATACCAACCACT CAAAAAG 

******:+:** i**^************:*:*********:*:*************** 


Mt: A AACCGATGGACT AT AAT AGT-T CT CCGCCAACAT ATGG AT C 

^0• ** *_(^**** ****CG*CC*ACGAAAAGT*A*************** *** 

fdq • _CA*CC*ACGAATAGT *A****************** 


Mt: AAAGACAAACTAT-CTTGCAAAGAAGCTGTCA 

Me: ******* ******CTACCA****************** 

Mg; *************CTGCCA****************** 

Figure 5. Comparison of nucleotide sequences of the variable region of the adhesive protein genes of 
Mytilis irossulus. M eduhs. and M galloprovincialis. Sequences of M cdulis and M galloprovincialis were 
according to Filpula el al (1990) and lnoue and Odo (1994), respectively. Asterisks indicate nucleotides 
identical to those of the M irossulus sequence. Hyphens indicate gaps inserted to align sequences. Ml, M 
Irossulus; Me, M eduhs; Mg. M galloprovincialis. 
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Figure 6. Amplification of the variable region of the adhesive protein 
gene of wild and cultured mussels. Amplified products were electropho- 
resed on 4.8% NuSieve GTG agarose gel (FMC). Lanes 1-3, wild mussels 
collected at Tromso; Lanes 4-6, mussels cultured in Brittany; Lanes 7- 
9. wild mussels collected at Hiura; Lanes 10-12, mussels cultured on 
the Mediterranean coast of France. M. molecular marker (pUCl9 DNA 
digested with I lap II). 

exhibited two fragments are presumed to be hybrids 
between the two species. 

Discussion 

Among the five species of the genus Afytilus. AL edulis, 
AL galtoprovincialis, and AL trossulus have been called 
the “AL edulis complex. " Since they are morphologically 
similar and shell shape is often influenced by local envi¬ 
ronment, it is difficult to identify these species by mor¬ 
phological characteristics. Recently, allozyme characters 
have been used to clarify the taxonomy of these species 
(Koehn et al, 1984; McDonald and Koehn, 1988; Varvio 
eta/., 1988; McDonald et al, 1991; Coustau et al, 1991; 
Viard et al., 1994). These characters are recognized as 
reliable markers, but data for multiple loci are required 
for accurate identification of all three species. Identifica¬ 
tion using mitochondrial DNA (mtDNA) sequences has 
also been described (Edwards and Skibinski. 1987; Blot 
et al.. 1990; Geller et al., 1993, 1994). Although such 
attempts were partially successful, it is still difficult to dif¬ 
ferentiate the three species unambiguously. In this study, 
we found that differences in a certain “variable region” 
of a sequence in the nonrepetilive domain of the foot 
protein 1 agree well with the taxonomic rank of species. 
It was also shown that the variations can be attributed to 
differences in the length of the fragments amplified by 
PCR. Thus the variable region may become an effective 
diagnostic marker. Because PCR requires only a small 
amount of DNA as a template, the method may be used 
for larvae or young individuals that are too small for 
analysis by other methods. In addition, the PCR system 
seems effective for the identification of hybrids within the 
AL edulis complex—we could detect hybrids that have 
the haplotypes of both Al. galloprovincialis and AL tros- 


sulus (Fig. 6). This system may become a powerful tool 
for studying the distribution and genetics of mussels, one 
of the most cosmopolitan of marine animals. 

We also tried to amplify the adhesive protein gene of 
AL coruscus, a mussel species native to Japan, but we 
were unsuccessful (data not shown). The adhesive protein 
gene sequence of Al. coruscus may be considerably dif¬ 
ferent from those of species in the AL edulis complex. It 
seems that other primers are required for analysis of mus¬ 
sels other than the Al. edulis complex. 

In contrast to the nonrepetitive domain, the length of 
the repetitive domain was highly variable, even in samples 
collected at the same site. This result is consistent with 
the fact that repeat patterns observed in three nucleotide 
sequences encoding foot protein 1 of AL edulis reported 
separately (Strausberg et al.. 1989; Filpula et al.. 1990; 
Laursen, 1992) differ from one another. Because the length 
of this region was partially overlapping among species, it 
seems inappropriate for use as a marker of species. Many 
individuals of the three species also have two or more 
fragments of different length. These results suggest that 
each species has enormous variability in the patterns of 
repetitive sequences. Considering that foot protein 1 is a 
key molecule for adhesion of mussels and thus is essential 
for their survival, its extensive diversity is intriguing. We 
suggest that the repetitive domain has been differentiated 
by a “shuffling” of the repeat pattern as well as by the 
common base substitution process (Inoue et al, 1995b). 
Increased diversity may have an important role in mussel 
survival or evolution. 
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